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Abstract 


We present in this paper a quantitative method for defining void size 
in large-scale structure based on percolation threshold density. Beginning 
with two-dimensional gravitational clustering simulations smoothed to the 
threshold of nonlinearity, we perform percolation analysis to determine the 
large scale structure. The resulting objective definition of voids has a natural 
scaling property, is topologically interesting, and can be applied immediately 
to redshift surveys. 


I. INTRODUCTION 

There has been increasing interest in large voids in the galaxy distribution. However, in 
large-scale theory and practice, there currently exist many different definitions of a void (for 
discussion see Sahni, Sathyaprakash &; Shandarin, 1993). Our purpose here is to present 
an objective and quantitative method for defining voids. We use an algorithm developed 
by Kauffmann (see Kauffmann & Melott 1992, hereafter KM) to locate voids. We use a 
percolation technique suggested by Shandarin (1983) in the form developed by Klypin & 
Shandarin (1993, hereafter KS) to determine the threshold density. We perform the analysis 
on two-dimensional N-body generated density distributions used in the work of Beacom et 

al (1991, hereafter BDMPS). 

The choice of the percolation threshold is motivated by the fact that it marks a change of 
topology. In the absence of percolation, the regions above the density threshold are isolated. 
On the contrary, in two-dimensional (hereafter 2-D) space, if the fraction above the threshold 
percolates, the fraction below the threshold forms isolated voids. 

We perform this work in two dimensions because coding and de-bugging are obviously 
simpler and visualization is easier. The agreement between the 2-D work of BDMPS and 
Melott & Shandarin’s 3-D work (1993, hereafter MS) is very impressive, suggesting 2-D 
is a good guide to what will happen in 3-D. The technique can be easily applied to two- 
dimensional galaxy samples, as well as to the microwave temperature fluctuation maps. 
After testing this method in 2-D we plan to study more realistic 3-D distributions. 
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II. NUMERICAL SIMULATIONS 

The numerical models are the ones used in the wnrlr nf rdmpc t1o2 , 

three-dimensional fi = 1 FRW^ u id verse T", vale . nt to a cross-section of a 

see BDMPS TheraT-? 3 ’ reS , P ? ctlvely ! ln three dimensions. For more detailed information 
The l m t US ‘sp^rnm 

Fbr each°type 

&s£tfr~’ *"■ " 1 “ 


D 2 (t) J kNL P(k)d 2 k = 1, 


universe «m i. "“t ° fgravitational instability, D(t) = a{t) for the fi = 1 

terms a(t l “*/* T ^ r ^ ^ ^ “ the absence oi pressure and radiation 
2 3 ’ 7 and k - 2* IT n °”' bneanty fre< l uenc y 15 c h°sen to be = 2 r k f where r = 

JSr Th - "t*" ■ “-i 


III. SMOOTHING 

We smooth each density distribution by Fourier convolution with the Gaussian window 

W(R) = e-*V»!. 

^ } - “-f- 


PERCOLATION 


sy, mg it as a tool to study topological properties of non-linear density fields. 
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n e ®! ence > we use percolation to study phase transitions also. In our case however it is 
not the system under study that is evolving, we have already evolved our simulations to 
the various points m time or epochs we wish to study. Rather, we take each of our ’snap 
shots of gravitational clustering and allow the ’threshold density’ to change. We choose 

th^ov^V^ 6 ’* 0 W j Cn j lty and mcrement or decrement this density unt/percolation in 
Tti ; J l e r d f n t Se or “ nd ® r - dei ? se Phase is reached using a square-lattice percolation algorithm. 
The percolation threshold is the point at which the transition is made from discontinuous 

threshold ° a globa cluster which spans the simulation box. The density at which this 
the smooth h 15 < :aJled the percolation density. This percolation density is applied to 

the smoothed density files such that anything equal to or above percolation density becomes 

simrd- nd t any hln ,? b t, OW , becon ? eS a ?• We caU this our percolated file. For the sake of 
simplicity, we will call regions above the threshold ’superclusters’ and those below it ’voids.’ 

thl„ A that , these ar ? I ! ot necessarily identical to the observationally based use of 

these words, and emphasize this by using single quotes. 

(i. In , tbe lde J* continuous 2-D world, percolation in the over-dense phase would mean 
e absence of percolation in the under-dense phase and vice-versa. However we study 

oercolate° n l” a ?' D and therefore « may be possible that neither phase 

DCTcolatfon thl °h d m - ° *T ake th ! P rocedure symmetric, we apply the void search at the 
♦ P h P ,lrrTr e !n° d u ln ,. t i he ° i T r ‘ den J Se p , hase and > measuring the ’supercluster’ sizes, we use 
the P ercolatl °n thresholdin the under-dense phase. This technique can be used to study 
n * volds , w j len the - v are isolated, since the over-dense phase percolates and vice 
and Shandl^(1992) a m ° re ^ dlscusslon of Percolation in 2-D systems in Dominik 


V. VOID SEARCH 

a W f perform a ’void’ search on the percolated file. The ’void’ search algorithm (KM) first 
finds the largest square of empty cells in the percolated file. Any empty cells along its four 
cannot r b% n added to tb is base void’, subject to the criterion that the length of the fill-in 
cannot be less than two-thirds the base ’void’ length. When the fill-in is complete, the next 

3, S ‘ 3 are 15 fo f T dand «*“?"« added t.° it- This process continues for smaller and 
smaller^ squares until all the voids m a particular simulation are found. We weight each 
void size by its predominance in the simulation. We characterize each simulation by the 
diameter of a circle of area equal to the weighted average ’void’ size y 


< D> = 


VdW d 

VW d 


where the summations are over the range of ’void’ sizes, d = 2,3,..., dMAX and 

N d d 2 


W d = 


d\iAX 


is the weight of ’voids’ of size d. N d is the number of Voids’ of size d. For purposes of 

soT&TtT , we . | erfor T d a supercluster’ search by percolating along under-dense regions 
so that the void search was actually finding the typical sizes of the ’superclusters ’ 8 


VI. RESULTS AND DISCUSSION 


We present the results in the form of the percentage area occupied by ’voids’ or ’super- 
clusters as a function of the diameter in figures 1-3, for n = 2,0, -2 respectively. Panel (a) 
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in figure 1 and 2 shows the distribution for the Gaussian fields with the corresponding power 
spectra. One can see that in Gaussian fields there is no statistical difference between the 
void and supercluster’ distributions, which was of course anticipated. It is worth stressing 
that the total area in ’voids’ and ’superclusters’ is only % 12% in the n = 2 model and 
% 31% in the 7i = 0 model. It reflects the fact that the ’void’ search algorithm can measure 
only large ’voids’ (see figure 4). One can see also that in the range D ~ 0.005, there are 
about 2.5 times more ’voids’ in the n = 0 model than in the n — 2 model and there are a 
lew larger voids as well. Both features are in agreement with the greater smoothness of 
the n = 0 model. 

The non-linear stages shown in panels b through g of figure 1 and 2 show the major 
difference between the ti = 2 and 7i — 0 models; in the n ~ 2 model ’superclusters’ occupy 
more area than ’voids’ and the largest structures are also ’superclusters,’ in the n = 0 model 
the opposite is true. The degenerate n = —2 model demonstrates that ’voids’ are always 
the dominant structures. 

The n = 2 model shows that both the ’voids’ and ’superclusters’ roughly scale with 
the scale of non-linearity A and the n = — 2 model does not, as can be seen in figure 1 
through 3 and especially in figure 5.. The mean diameter of ’superclusters’ is slightly but 
significantly greater than that of ’voids’ in the ti = 2 model, and both are about 2 times 
smaller than Ajvx, shown as a dotted line in figure 5. 

In the n = 0 model ’voids’ are a little larger than ’superclusters’ and the difference with 
Ajvx is somewhat less, about 1.8 times. The models also do not display quite as good scaling 
with X NL as the n = 2 model. 6 6 


As previously mentioned in the section on the void search routine, we don’t necessarily 
include all of the area below the percolation threshold in ’voids.’ The ’void’ search routine 
was written to approximate circular ’voids.’ This is most obvious in figure 4, where the 
early epochs show a combined ’void’/cluster area of considerably less than 100%. In the 
later epochs, the voids’ become larger and more circular; less ’void’ area is lost due to the 
two-thirds requirement during fill-in (see description of void search). The careful observer 
will also notice that a few of the data points are slightly above 100%, most noticably the 
last two data points on the n = —2 (dashed) line. This is simply because we measure 
voids at the percolation threshold of the over-dense regions and ’superclusters’ at that of 
the under-dense regions, and these percolation thresholds are not identical. 

The advantages of our approach are both theoretical and practical. The smoothing 
we use removes discreteness effects and our choice of smoothing length is insensitive to 
statistical fluctuations. The use of percolation objectively defines a connected structure, and 
its combination with an objective measure of void size remedies its major shortcoming as a 
large-scale structure statistic, i.e., it does not give rise to a physical lengthscale associated 
with percolating structures. As shown in earlier work, the results of percolation studies can 
be connected with the initial conditions that generated the structures. 
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IX. FIGURE CAPTIONS 

Figure la) Percentage area occupied by Voids’ (solid lines) and ’superclusters’ (dashed 
lines) of size L (units of box size) for a purely Gaussian density distribution with k c = 256k f 
m the n = 2 model. 1 

Figure lb) Percentage area occupied by ’voids’ (solid lines) and ’superclusters’ (dashed 
lines) of size L (units of box size) for the epoch k NL = 128 in the n = 2 model. 

Figure lc to lh) Same as Fig. lb but k^i as specified on plot. 

Figure 2) Same as Fig. 1 for n = 0 model 

^9 S i 1 n e fi 3a ^? ame t S Fig ‘ 1 for n = “ 2 model except that only stages k NL = 
o4, 32, 16, 8, and4 are shown. 

/ areas occu pied by both ’voids’ and ’superclusters’ versus X NL for n — 2 

(sohd line), n = 0 (dotted Une) and n = - 2 (dashed Une) models. 

Figure 5a) Evolution of mean ’void’ size (solid line) and ’supercluster’ size (dashed line) 
versus a^l in the n — 2 model, d = \ NL (dotted line) is included for comparison. 

Figure 5b) Same as Fig. 5a but in the n - 0 model. 

Figure 5c) Same as Fig. 5a but in the n = -2 model. 
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